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I.  INTRODUCTION 


Continuous-  wave  (cw)  chemical  lasers  are  inhomogeneously  broadened. 
Under  lasing  conditions,  the  gain  medium  is  in  a  state  of  both  rotational 
and  translational  nonequllibrlura. 

Early  cw  chemical  laser  models  Incorporated  the  assumption  of  rota¬ 
tional  and  translational  equilibrium  and  gave  reasonable  estimates  of  net 
laser  output  power  and  scaling  laws  (e.g..  Refs.  1  and  2).  The  later 
inclusion  of  rotational  nonequilibrium  in  these  models  provided  estimates 
of  the  output  from  various  rotational  energy  levels  (e.g.,  Refs.  3-5).  The 
Inclusion  of  translational  nonequilibrium  (i.e.,  inhomogeneous  broadening 
with  hole  burning)  provided  estimates  of  the  output  from  each  of  the  multi¬ 
ple  longitudinal  modes  in  large  scale  cw  chemical  lasers  (e.g.,  Refs.  6-8). 

The  simultaneous  consideration  of  both  rotational  and  translational 
nonequilibrium  in  cw  chemical  lasers  has  received  less  attention.  A  gen¬ 
eral  analytic  model^  and  a  computer  code^*^’^^  have  been  developed  with 

limited  results  presented.  More  recently,  an  analytic  model  was  devel- 

1  “y 

oped  that  provides  solutions  for  both  amplifiers  and  Fabry-Perot  (F-P) 
oscillators . 

The  purpose  of  the  present  report  is  to  consider  both  rotational  and 
translational  nonequilibrium  and  provide  a  multiline  (multiple  rotational 
energy  levels)  and  multimode  (multiple  longitudinal  modes)  solution  for  a 
large  scale  cw  chemical  laser  with  an  F-P  resonator.  It  is  assumed  that 
Av  <<Av  <<Av  where  Av  ,  Av,  ,  and  Av .  are  frequencies  that  characterize 
longitudinal  mode  separation,  homogeneous  width,  and  Doppler  width,  respec¬ 
tively.  The  present  report  may  be  viewed  as  an  application  of  Ref.  12  that 
generalizes  Ref.  7  to  include  multiple  rotational  energy  levels.  The  nota¬ 
tion  is  the  same  as  that  in  Ref.  12.  Symbols  and  nondimensional  variables 
are  summarized  in  Appendixes  A  and  B,  respectively.  The  corresponding  case 
of  rotational  nonequilibrium  and  translational  equilibrium  is  discussed  in 
Appendix  C. 
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II.  THEORY 


The  present  model  is  briefly  described.  Laser  equations  are  for¬ 
mulated,  and  solutions  are  then  presented  for  a  cw  chemical  laser  oscilla¬ 
tor  with  laminar  mixing. 

A.  FLOW  MODEL 

A  cw  chemical  laser  with  an  F-P  resonator  is  illustrated  in  Fig.  la. 
The  present  simplified  mixing  model  is  illustrated  in  Fig.  lb.  The 
reactants  are  ar-  umed  to  be  premixed  but  do  not  react  until  a  flame  sheet, 
yf(x),  is  reached.  The  flame  sheet  shape  is  specified,  a  priori,  from 
diffusion  theory.  The  streamwise  station  where  the  flame  sheet  reaches  the 
channel  center  line  is  denoted  x^  and  characterizes  the  diffusion  rate. 
The  width  per  semichannel  and  the  number  of  semichannels  are  denoted  w  and 
n^^,  respectively.  Laser  radiation  is  in  the  ±y  direction. 

B.  GAIN  AND  RESONATOR  MODEL 

In  order  to  simplify  the  mathematical  development,  we  assume  a  "0" 
type  laser  transition 

v+l,J-*-v,J  (la) 

rather  than  the  "P"  type  laser  transition 

v+l,J-l-^v,J  (lb) 

appropriate  for  cw  chemical  lasers.  Here,  v,  J  denotes  vibrational  and 
rotational  energy  level,  respectively.  This  Q  type  laser  transition 
approximation  is  used  in  Ref.  11  and  is  consistent  with  other  simplifying 
assumptions  used  in  the  present  model.  In  addition,  we  consider  a  two- 
vibrational  energy  level  model  and  denote  the  upper  and  lower  levels  by 
subscripts  2  and  1,  respectively.  Expressions  for  gain  are  given  by  Eqs. 
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(B-5a)-(B-5d) .  The  quantity  n^(J,v),  In  these  equations,  denotes  the  num¬ 
ber  density  of  particles  in  the  v,  J  level  that  are  resonant  with  radiation 
in  the  frequency  range  v  to  v  +  dv.  (In  the  present  study,  as  in  the  pre- 

C.  10 

vious  studies,  resonant  frequency  is  used  instead  of  particle  thermal 

velocity  in  order  to  evaluate  the  Interaction  of  the  radiation  field  with 
the  Inhomogeneously  broadened  gain  medium.  For  radiation  in  the  ±  y  direc¬ 
tion,  the  resonant  frequency  v  is  related  to  the  particle  thermal  velocity 
Vy  by  the  Doppler  relation  “  1  =  +Vy/c.  Similarly,  n^(J)  and  n^ 

denote  the  net  number  density  of  particles  in  the  energy  levels  v,  J,  and 
V,  respectively.  For  the  case  of  translational  and  rotational  equilibrium, 
the  quantities,  n^(J,v),  n^(J),  and  n^  are  related  by  Eqs.  (B-4a)- 
(B-4d).  Note  from  Eq.  (B-Ad)  that  f,  denotes  the  fraction  of  particles  in 

*J 

rotational  level  J  for  the  case  of  rotational  equilibrium.  The  maximum 

value  of  f ,  is  denoted  f,  and  is  a  function  of  On  =  T^/T  (Table  I).  It 

J  J,m  ^  X 

is  convenient  to  let  f  =  f,  ,  where  f  is  the  reference  value  of  f,  used 

r  J,m’  r  J 

in  the  normalizations  in  Eqs.  (B-3a)-(B-3d). 

The  two  end  mirrors  in  the  F-P  resonator  are  each  assumed  to  have  the 

same  reflectivity,  R^.  Line  center  laser  frequency  is  denoted  [e.g., 

Eq.  (B-Ac)],  and  the  dependence  on  J  is  ignored.  The  mirror  separation  L 

is  adjusted  so  that  the  longitudinal  mode  frequencies  v  .  are  symmetric 

d  »  J 

about  V  .  The  longitudinal  mode  frequencies  are  then 
o 

V  =  V  -t-  j  Av  i  =  0,  ±1 ,  +2 , . . .  (2a) 

J,j  o  c 

or 

''j,j  “  ^o  "  T^^'^c  j  "  1.  2,  3,...  (2b) 

=  ^o  (-1  +  .1  =  -1.  -2.  -3,... 

where  j  and  Av^  denote  longitudinal  mode  number  and  longitudinal  separa¬ 
tion,  respectively.  Let  and  denote  local  radiation  Intensity 

d  d  » J 

in  the  +y  and  -y  directions  (Fig.  lb).  Due  to  the  symmetry  of  the  longi- 
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Table  I .  Maximum  value  of  fj  and  corresponding 
value  of  J  for  fixed 

H 


Jm 

f, 

J  ,m 

1 .  1 0 

2 

0.265A 

1/20 

3 

0. 1889 

1  /AO 

A 

0. 1353 

1/60 

5 

0.1106 

1/80 

6 

0.0957 

1/100 

7 

0.085A 

^Eq.  (B-Ad). 
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tudinrJ  modes  about  v  the  interaction  of  I 

o  ’ 

mecl^um  can  be  evaluated  by  consideration 

as  discussed  In  Ref.  6. 
intensity  with  frequency  is  illustrated  in  Fig 
highest  frequency  at  which  lasing  occurs. 


^  ^  and  with  the  gain 

»  J  ^  J 

of  the  net  local  intensity 

A  typical  variation  of  local 

.  2.  Here,  v  denotes  the 
J  ,  t 


C.  LASER  EQUATIONS 

Consider  the  case 

R  »  1;  R  »  1 
t  ’  r 


^^h  ^  p(Torr) 

Au,  100 
d 


Av 


50 


Av,  L(ra)p(Torr) 


«1 


1  - 


Av^ 
Av , 


«  1 


(3a) 


(3b) 


(3c) 


(3d) 


where  L(m)  is  resonator  mirror  separation  in  meters,  p(Torr)  is  gain  region 

pressure  in  Torr,  and  where  R^  =  and  denote  the  ratio 

of  the  translational  and  rotational  relaxation  rate  to  the  rate  of  colli- 

1 2 

sional  deactivation.  In  general,  R^'^^r  ~  0(1)  cw  chemical  lasers. 

Equations  (3b)  and  (3c)  imply  multiple  longitudinal  modes.  These  equations 
apply  for  large  scale  chemical  lasers^’^  (e.g.,  p  =  10  Torr,  L  =  50  ra)  and 
are  Illustrated  in  Fig.  3.  As  a  consequence  of  Eqs.  (3a)  and  (3d),  the 

equations  that  define  the  performance  of  a  cw  chemical  laser  with  an  F-P 
12 

resonator  are 


Jv 


O  ^ 

-X"N 


1  +  0  {(1  -  (R^/R^)l(Av^/Av^)} 


1 1 


(4a) 


Fig.  2 


l’  -  v  _■  t'-  i'  -  2  f'j  • 


Local  Longitudinal  Mode  Intensity  in  F-P  Resonator  for 
Case  Where  v,  ,  Is  Symmetric  About  v  . 

-  (-^)  .  T  (-) 


Note : 


I 


J  ,  V  ^  J  »  V 


I 


J.J 


[Eq.  {2a)]. 
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Fig.  3.  Characteristic  Frequencies 
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(4b) 


dN'/di;  =  dN  /d^  -  N_  -  N  -  (SG^/ir)  III,. 

•L  ^  c  .  J  ,  .1 


Gj(V) 


l(v  -  u') 


dv 

Av, 


(4c) 


where  is  the  threshold  gain,  and  other  variables  are  defined  in 
Appendixes  A  and  B.  These  equations  can  be  evaluated  without  consideration 
of  Nj,  which  is  found  from 


r^Jln2^1/2  ^^h  ^t  V  X 

-  ■  aT:  r~  7^  4  1 

d 


r  fjN 


J.j 


(4d) 


Here,  N,  ,  N  ,  N  ,  N  ,  and  I  .  are  normalized  forms  of  the  variables 

J  »  J  Jl  d  ^  J 

n^CJ.v)  -  nj(J,v),  n2(d)  -  n^CJ),  1x2  -  Hj,  n2  +  nj ,  and  respec¬ 

tively.  The  quantity  I  .  has  been  normalized  by  a  saturation  intensity, 

_  d  »  j 

G^k^^/(2a^),  which  characterizes  hole  burning  (translational  nonequil- 
ibrium)  effects.  I^Then  I,  ,  is  small,  the  laser  medium  is  in  translational 
and  rotational  equilibrium.  Equations  (la)  and  (Id)  then  indicate 


N 


Jv 


f 

f^e  N 


=  1 


5) 


When  I  is  not  small,  a  nonequilibrium  solution  of  Eqs.  (4a)-(4d)  ls 
J » .1 

required . 


Equation  (3c)  indicates  that  there  are  many  longitudinal  modes  within 

the  homogeneous  width  Av.  ,  In  order  to  facilitate  the  summation  of  I,  . 

h  J » j 

L(v  -  V  )  in  Eq.  (4a),  it  is  convenient  to  replace  I,  .  by  a  continuous 
J  ,  j  ^  .  I 

distribution  I,(v)  in  the  interval  from  v,  .  -  Av  /2  to  v,  .  +  Av  /2.  Thus 
J  J,j  c  J,j  c 


V,  .  +  Av  /2 

/  I,(v)dv  =  I.(v  .)Av 

J  J  J ,  1  c 

V  -  Av  /2 

c 


(6) 
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It  follows  that 


=  (Tr/2)Av^Ij(v) 


(7a) 


■  I  'j.j  'i 


(7b) 


J  ^ 


(7c) 


Let  denote  the  largest  lasing  frequency  for  a  given  J.  For  the 
present  case  of  an  F-P  resonator  with  closely  packed  longitudinal  modes 
(etg.,  Av  <<  Av  ) ,  it  is  expected  that  N  will  depart  only  slightly  from 

C  n  J  V 

the  threshold  value  for  frequencies  in  the  range  jxj  <  Xj^  (Fig.  4). 

The  resulting  particle  density  distribution  can  then  be  approximated  by 


1 

„2  2 

1*1  <  *jf 

(8a) 

X,.  -  X 
=  e  Jf 

1*1  >  *jf 

(8b) 

which  assumes  that  is  continuous  at  X  =  X^^  and  has  an  equilibrium 


distribution  for  |X|  >  Substitution  of  Eqs.  (8a)-(8b)  into  Eqs.  (4a)- 

(4d)  and  (7a)-(7c)  yields 


(2/1.  )G^  -  Njv, 


(9a) 


2  2 

(Tr/2)Av^  *  e  -1 


x|  <  X  (9b) 


(TT£n2)^^^  ®  ^jf  ■  (2/r^^^)Xj^ 


(9c) 
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16 


2 

*Jf  “  ^  f(2/n)G 


■] 


where  it  is  assumed,  in  Eqs*  (9b)  and  (9d),  that  ■  0  at  |xj 
Equation  (9b)  indicates  an  exponential  variation  of  r^(v)  with  X. 

s/ 

Output  power  per  unit  volume  is  found  from  [e.g.,  Eq.  (B-6c)] 


(9d) 


Jf ' 


dPj(v)/dz;  -  [SG^/(2Tt)]tj(v) 

(10a) 

OD 

dPj/d;  =  [SG^/(2it)]Ij  =  /  Pj(v)dv 

-»oo 

(10b) 

dP/d^  =  y  dPr/d? 

J 

(10c) 

Integration  of  Eqs.  (lOa)-(lOc)  to  the  end  of  the  lasing  region  5  provides 
the  corresponding  net  output  values  [P^(v)]  ,  P,  ,  and  P  . 

J  0  J  1 6  0 


D.  LAMINAR  MIXING 

fi  12 

For  the  case  of  laminar  mixing  » 

N^  = 

Eq.  (4b)  may  be  expressed  in  the  form 

Zp(dN"/dz)  -  1  -  2Zj^zN~  -  2z^  -  2z4,^  I  (iTin2)^'^^( AVj^/Avj)Ij  (11) 

3 

where 


,1/2.  Z  .  c  1/2 

^  ’  D 

(12a) 

(12b) 
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Integration  of  Eq.  (11)  for  the  case  of  zero  power  yields 


where 


1 

is  the  Dawson  integral, 
where  lasing  is  initiated. 


ZjjN~  “  2D(z)  -  z 


(12c) 


D(z) 


(12d) 


Let  subscript  i  denote  conditions  at  the  station 
The  boundary  condition  at  this  station  is 


z  N  =  (2/tt)z  G  /f  =  2D(z  )  -  z 
Di  DcJ,m  1  1 


(13) 


where  f j  denotes  the  maximum  value  of  fj.  Corresponding  values  of  z^, 

z^^N.  ,  and  (2/'n)z„G  /f,  ai"®  listed  in  Table  II  and  are  used  interchange- 
D  1  D  c  J  ,ra 

ably.  It  is  convenient  to  choose  f  =  f.  ,  so  that  fr  _  =  1.  However, 

r  j  ‘j 

for  generality,  will  be  retained  in  subsequent  equations.  The  inte¬ 

gration  of  Eq.  (11),  subject  to  the  boundary  conditions  given  by  Eq.  (13), 
is  now  discussed  for  equilibrium  and  nonequilibrium  cases. 


1 .  EOUILIBRIITM  (R^  ->•  ») 

In  the  limit  of  translational  and  rotational  equilibrium 
(Rj  R^  ■*■  “),  ->•  0  and  lasing  occurs  on  line  center  at  the  single 

rotational  level  corresponding  to  the  maximum  value  of  fj.  Thus,  in  the 
region  Zj  ^  z  ^ 


Also 


z„N"  =  z„n7  =  (2/tt)z„G  /f, 

D  Di  'DcJ,m 


(Av./Av,)I  =  (1  -  2zz_n7  -  2z^)/z 
t  n  a  u  1 


(14a) 

(14b) 


=  2D(z^)  - 


^D^i 


(14c) 
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TABLE  II.  Relation  between  station  where  lasing  is  initiated  z^, 

threshold  number  density  N",  and  threshold  gain  G. 

X  c 

for  laminar  mixing.^ 


0.001 

0.0010 

0.005 

0.0050 

0.010 

0.0100 

0.050 

0.0498 

0.100 

0.0987 

0.200 

0. 1895 

0.300 

0.2653 

0.400 

0.3199 

0.500 

0.3489 

0.552 

0.3528^^ 

®(Eq.  (13). 

^Maximum  zero  power  value. 
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3i  e 

z. 


(l^d) 

( 14e) 


(15) 


^'o^'e  ■  [*  ■  ■  (2/3)z^] 

In  the  further  limit  -►  0, 

-►  0  ;  >  2“^^^  ;  z^^P^  >  (18)"^^^ 

which  corresponds  to  saturated  laser  operation. 


2 .  NONEQUILIBRIUM  (R^  »  1 ,  »  1) 

Under  conditions  of  translational  and  rotational  nonequilibrium ,  a 
numerical  integration  of  Eq.  (11)  is  generally  required. 

For  z  near  z^ ,  the  last  term  in  Eq.  (11)  is  negligible,  and  the  zero 
order  solution  is  applicable.  Thus  the  leading  terms  of  the  solution,  near 
z  =  z^,  are 

1  -  4z .D(z. ) 


li-  1  =  r 

Nj  ^  2D(z^)  -  z. 


■]  (z  -  z^) 


=  [(n”/n’)  -  1] 


1/2 


(16a) 

(16b) 


Note  that  (N  /N\  )  -  1  has  a  linear  dependence  on  z  -  z^ ,  while  Xjf  varies 
as  (z  -  z^)^^^. 

In  order  to  integrate  Eq.  (11),  it  is  necessary  to  specify 
6j^,  f^,  R^ ,  and  z^  (or  Zj^N^).  The  following  quantities  can  then  be  found; 


ZdN-  =  2D(z^)  -  z^ 


=  z_,N,  f-  f 
t  D  1  J,ra  r  t 


fj  (2.1  +  1)  expt-J(J  +  1)9^^] 
fT7  "  (2J  +  1)  exp[-J  (J  +  l)9p] 


(17a) 

(17b) 

(17c) 


m  m 
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The  integration  of  Eq .  (11)  proceeds  by  the  evaluation  of 


Xjj  .  (17d) 

(TT£n2)^'^^(Av^/Au^)Ij  -  e  erf  -  (2/ir^^^)Xj^  (I7e) 

at  successive  streamwise  stations.  The  range  of  J,  at  each  strearawise  sta¬ 
tion,  is  limited  by  the  requirement  that  Xj£  be  real  in  Eq.  (I7d).  Net 
output  is  found  from  Eqs.  (lOa)-(lOc).  Thus, 


o 


2  2 
z  X^  -  X^ 

^  (e  -  Dzdz  (18a) 


'/j,e  -  =<Jf  --[72 

O  TT 


(18b) 


Ve  =  1  C"') 

where  the  integrands  in  Eqs.  (18a)  and  (18b)  are  either  nonnegative  or 
zero.  (That  is,  the  effective  integration  interval  in  Eqs.  (18a)  and  (18b) 
corresponds  to  values  of  z  for  which  the  integrand  is  nonnegative.) 
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III.  NUMERICAL  RESULTS  AND  DISCUSSION 
Solutions  of  Eqs.  (17)  and  (18)  have  been  obtained  for  cases  where 


=  0.001  -  0.5  (19a) 
Qj^  =  1/80  (J„  =  6,  fj^^  =  0.0957)  (19b) 
^r  =  fj.m  =  0.0957  (fj^^  =  1)  (19c) 
Rj.  =  100  (19d) 
R^  »  1  (19e) 


Results,  which  are  presented  in  Figs.  5  to  8,  are  discussed  herein. 

Figure  5a  indicates  the  variation  of  inversion  number  density  z^N 

with  streamwise  distance  in  the  lasing  region  z^  <  z  <  for  a  variety  of 

threshold  gains.  (Recall  from  Table  II  that  and  (2/ii)2pG^  are  used 

interchangeably).  The  inversion  number  density  z^^N  increases  above  the 

threshold  value  z^N.  for  z.  <  z  <  z  .  In  the  limit  R  *  <*>,  z_N  =  z,.N.  in 
Di  1  e  t’DDi 

this  region. 

The  variation  of  local  laser  intensity  Ij  with  streamwise  distance  is 
given  in  Fig.  5b  for  the  case  z^  =  0.1.  In  accord  with  Eq.  (19b),  maximum 
intensity  is  obtained  for  J  =  6.  Note  that  for  the  present  case  of  an  F-P 
resonator,  the  peak  intensity  for  each  J  level  occurs  at  the  same  stream- 
wise  station.  The  streamwise  extent  of  the  lasing  region  decreases  with 
departure  of  J  from  J  =  6. 

The  number  of  longitudinal  modes  at  a  given  streamwise  station  is 
approximately  equal  to  2  Xj^/Av^.  The  variation  of  Xj^  and,  therefore,  the 
variation  of  the  number  of  longitudinal  modes  with  streamwise  distance,  is 
given  in  Fig.  5c  for  =  0.1.  It  is  clear,  from  Figs.  5b  and  5c,  and  from 
physical  considerations,  that  Xjf  is  a  maximum  where  Ij  is  a  maximum  and 
that  the  streamwise  extent  of  Xj^  is  the  same  as  that  for  Ij. 
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STREAMWISE  DISTANCE  2 


Variation  of  Inversion  Number  Density,  Laser  Intensity,  and  Longi 
tudinal  Mode  Frequency  Range  for  F-P  Resonator  with  Operating  Con 
ditions  Given  by  Eq.  (19).  (a)  Inversion  number  density  in  lasing 
region,  2.  <  2  <  z  . 


LOCAL  LASER  INIENSIfY.  I, 


LONGIIIJDINAI  MnOf  FREQUENCY  RANGE 


.  =  6 


5 


=  '  8C 
R  =  ’OC 


(a) 


-!  e  e  'c' 


Net  Output  Power  from  Each  Rotational  Level  Versus  (a)  Threshold 
Gain  Parameter  ZpNT  and  (b)  Rotational  Level  J  for  Operating 
Conditions  Given^by  Eq.  (19).  Note  that  ordinate  may  be  viewed  as 
fraction  of  net  output  power  from  a  saturated  laser. 
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The  vailatlon  of  longitudinal  mode  power  (P^  )  =  Av  [P  (v)]  with 

frequency  is  given  in  Figs.  6a  and  6b  for  =  0.001  and  0.01,  respective¬ 
ly.  These  curves  are  approximately  exponential  of  the  form  exp(-X^)  in 
accord  with  Eq.  (9b).  The  departure  from  the  exponential  form  is  due  to 
the  variation  of  Xjf  with  z  in  Eq.  (18a). 

Net  output  power  per  rotational  level  P.  ^  is  given  in  Figs.  7a  and 

J  ,e 

7b.  Maximum  power  is  obtained  on  J  =  6  as  expected.  The  variation  of  P.  „ 

J  ,e 

with  threshold  gain  is  given  in  Fig.  7a.  Note  that  the  power  Pt 

o 

increases  with  increase  in  threshold  gain,  reaches  a  maximum  at  approxi¬ 
mately  (2/tt)  =  0.05,  and  then  decreases  with  further  Increase  in 

threshold  gain.  The  reason  for  this  behavior  can  be  deduced  from  Figs.  7 

and  8.  Thus,  it  is  seen  from  Figs.  7b  and  8  that  the  number  of  active 

rotational  levels  J  decreases,  while  net  laser  power  remains  fairly  con¬ 
stant,  with  Increase  of  threshold  gain  in  the  region  where  the  threshold 

gain  is  small.  Hence,  for  the  present  cases,  the  output  on  lines 

3  <  J  <  9  increases,  while  the  output  on  lines  J  <  3  and  J  >  9  decreases 
with  an  initial  increase  in  threshold  gain,  as  indicated  in  Fig.  7a. 

When  (2/tt)  =  0.05,  the  output  at  each  J  level,  as  well  as  the  number 

of  active  J  levels,  both  decrease  with  further  increase  in  (2/ti )ZjjGj,. 

Figure  8  indicates  the  effect  of  threshold  gain  and  translational 

relaxation  rate  on  net  output  power.  When  (2/TT)ZpG^  <  0.05,  the  results 
for  the  realistic  value  =  100  agree  within  about  5%  with  those  for 

translational  equilibrium  (R^^  -►“').  Hence  equilibrium  solutions  tend  to 

provide  realistic  estimates  for  net  output  power  from  cw  chemical  lasers 
operating  under  saturated  conditions. 

The  corresponding  case  of  rotational  nonequilibrium  and  translational 

equilibrium  is  discussed  in  Appendix  C.  The  latter  is  a  common  assump¬ 
tion.  In  this  case,  there  Is  a  single  longitudinal  mode,  at  line 

center,  for  each  rotational  energy  level  with  sufficient  gain  to  lase.  The 

output  from  each  rotational  level,  and  the  net  output  power,  are  Indicated 

in  Fig.  9  and  Table  III,  respectively.  The  output  is  Intermediate  between 
the  case  of  rotational  and  translational  equilibrium  and  the  case  of  rota- 
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NFT  ourpul  PfR  ROTAIIONAL  LEVEL  (18)  / 


THRESHOLD  GAIN  (2'7r)  2p  G, 


Fig.  9.  Net  Output  Power  from  Each  Rotational  Level  for  Case  of 
Translational  Equilibrium  and  Rotational  Nonequilibrium 
[Eq.  (C-3)] 
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TABLE  III.  Net  output  power  for  cw  chemical  laser  with  laminar 

mixing,  9  =  1/80,  and  (a)  translational  and  rotational 

R 

Equilibrium,  (b)  translational  equilibrium  and  rotational 
nonequilibrium,  or  (c)  translational  and  rotational 


nonequilibrium  with 

Av  «  Av  « 
c  h 

Av ,. 
d 

Net 

Output  Power 

,  (18) 

1/2^  p 

zi 

(a)  Rj.  » 

(b)  R^  = 

100 

(c)  Rj.  »  1 

Rj.  -*■  » 

Rt  " 

OD 

R^  =  100 

0.001 

0.997 

0.986 

0.984 

0.005 

0.984 

0.871 

0.965 

0.010 

0.968 

0.953 

0.945 

0.050 

0.843 

0.822 

0.800 

0.100 

0.694 

0.671 

0.639 

0.200 

0.427 

0.405 

0.363 

0.300 

0.218 

0.  199 

0.156 

0.400 

0.077 

0.063 

0.033 

0.500 

0.009 

0.004 

0.000 

0.552 

0.000 

0.000 

0.000 
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tional  and  translational  nonequtlibrium.  The  differences  become  more  pro 
nounced  as  the  degree  of  saturation  is  decreased. 
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IV.  CONCLUDING  REMARKS 


The  present  model  provides  an  improvement  over  previous  models  (e.g., 
Refs.  3  to  8)  because  multiple  longitudinal  modes,  as  well  as  multiple 
rotational  laser  levels,  are  included.  Although  the  present  development 
assumes  an  F-P  resonator,  it  Is  expected  that  the  results  characterize  the 
performance  of  cw  chemical  lasers  with  more  realistic  resonator  configura¬ 
tions  . 

A  comparison  between  the  present  model  and  a  detailed  numerical  code 

calculation,  ^  ^  presented  in  Ref.  7,  suggests  that  Eq.  (3c)  is  overly 

restrictive  and  that  the  present  model  provides  reasonable  results  for 

Av  /Av,  =  0(1)  (e.g.,  p  =  10  Torr,  L  =  10  m) . 
c  h 
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APPENDIX  A.  PARTIAL  LIST  OF  SYMBOLS 


c 

DO 


G  (v) 
gj(v) 

Ij(v) 

J 

J 

*^cd’^tr'*^rr 

L(v-v‘ ) 
Nj;.Nj,n- 


n  (J,v),n  (J),n 

y'  ’  y'  ^  y 


Pj(v),Pj,P 


Rt  *Rr 
S 

T 


speed  of  light  in  vacuum 
Dawson  integral,  Eq.  {12d) 

fraction  of  particles  in  rotational  energy  level,  Eq. 
(B-i<d) 

reference  value  of  f,,  Eq.  (B-4d) 

J  r 

normalized  threshold  gain,  Eq.  (B-5d) 
normalized  gain,  Eq.  (B-5b) 
gain,  Eq.  (B-5a) 

intensity  for  longitudinal  mode  J,J 
net  intensities 
intensity  per  unit  v,  Eq.  (6) 
nondimensional  intensity,  Eq.  (B-6a) 
rotational  energy  level 
longitudinal  mode  number 

deactivation,  translational  and  rotational  relaxation 
rates.  Ref.  12 

Lorentzian  distribution,  Eq.  (B-5a) 

normalized  inversion  number  densities,  Eq.  (B-3) 

normalized  total  number  of  lasing  species,  Eq.  (B-3) 

particle  number  densities,  Eq.  (B-3) 

total  number  of  lasing  species,  n.|  +  02 

reference  value  for  1)2 

normalized  output  power  released  up  to  station  c,  Eqs. 
(10)  and  (B-6c) 

net  output  power  released  between 
stations  and  EQS.  (10)  and  (B-6c) 

reciprocal  Doppler  width,  Eq.  (B-4d) 

Fabry-Perot  resonator  mirror  reflectivity 

collisional  rate  ratios, 

parameter,  Eq.  (B-6c) 

temperature 
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APPENDIX  A.  PARTIAL  LIST  OF  SYMBOLS  (continued) 


X 

y 

yp 

z 

2d 

“^D 

X 


J,j 

V 

O 

Av 


characteristic  rotational  energy  temperature,  which 
equals  30.16K  for  HF,  Eq.  (B-4d) 

streamwise  velocity,  Fig.  la 

vibrational  energy  level 

thermal  velocity  in  y  direction.  Fig.  la 

normalized  frequency,  Eq.  (B-Mc) 

value  of  X  corresponding  to  largest  active  longitudinal 
mode  frequency  v,  Figs.  2  and  4b 

streaimwise  distance.  Fig.  la 

characteristic  diffusion  distance.  Fig.  la 

transverse  distance.  Fig.  la 

flame  sheet  ordinate.  Fig.  1b 


energy  per  mole  of  photons 

normalized  streamwise  distance,  x'u 

normalized  diffusion  distance,  Xj^/u 

characteristic  rotational  temperature  parameter,  Tj^.'T 

wavelength 

frequency 

longitudinal  mode  frequency 
line  center  frequency 

Doppler  width  [full-width,  half-maximum  (FWHM)] 

homogeneous  width  (FWHM) 

longitudinal  mode  separation,  c/2L 

cross  section  for  stimulated  emission,  Eq.  (B-5a) 

parameters,  Eqs.  (12b)  and  (C-lf) 
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APPENDIX  A.  PARTIAL  LIST  OF  SYMBOLS  (continued) 

Subscripts 
e 
i 
J 

j 
m 
r 

V 

Superscripts 

(+),(-)  radiation  in  +y,  -y  directions 

difference  of  number  densities,  Eq .  (B-3) 


end  of  lasing  region 
start  of  lasing  region 
rotational  level 
longitudinal  mode 
maximum  value 

reference  value  or  rotational  relaxation  value 
vibrational  level 
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APPENDIX  B.  NONDIMENSIONAL  VARIABLES 


Geometry 


y/w  »  (x/xjj)^/2  Laminar  Flame  Sheet 


5  x/u;  x^/u 


D  "cd  "D' 


Number  Density 


[n^CJ.v)  -  n^(J,v)]y^ 


J  V 


n  w  f  p 
r  r  *^0 


N,  = 


[n^CJ)  -  n^(J)]y^ 


n  w  f 
r  r 


n"  =  (n^  -  nj^)y^/(n^w) 


h'l  =  (>^2  “  n^yf/(nj.w) 


Equl librluTTi  Number  Density  Distributions 


_  _y2 

n^(J)  '  Po  ® 


Translational  Equilibrium 


P 


o 


r4£n2 


1/2 

)  /Av^ 


X  =  (4jln2)^^^  (v  -  V  )/ Av. 

o  a 


(B-1) 

(B-2) 

(B-3a) 

(B-3b) 

(B-3c) 

(B-3d) 

(B-4a) 

(B-4b) 


(B-4c) 


APPENDIX  B.  NONDIMENS lONAL  VARIABLES  (continued) 


Gain 


n^(J)  _  (2J  +  l)exp[-J(J  +  l)0p] 

"  l.(2J  +  l)exp[-J(J  +  1)6_]  (B-4d) 

=  Rotational  Equilibrium 


gj(v')  =  j  -  nj(J,v)]  L(v  -  v')dv  Q-Branch  (B-5a) 


L(X  -  X’)  =  {l  +  4[(v  -  v’)/Av^]^} 


-1 


Gj(v')  =  y^/(n^w  f^p^  Av^ 


=  /  n“^  L(v  -  v')dv/(Av^) 


(B-5b) 


g  =  (-Jin  R  )/(yrn  )  Threshold  Gain  (B-Sc") 

c  m  r  sc 


G 

c 


gcXf/CnrW  Av^  5^) 


(B-5d) 


H2 


APPENDIX  B.  NONDIMENS lONAL  VARIABLES  (continued) 


I 


Intensity  and  Power 


(B-6a) 


I  =  I  Ij  =  II  Ij  . 

J  J  3 


(B-6b) 


S  G  S 

P - P - /  Id; 

n  w  ue^  2it 
r  f  C, 


Av, 


(B-6c) 


=  Tip  f  R^ 

o  r  t  n 


Rates 


(B-7) 


M3 


APPENDIX  C 

ROTATIONAL  NONEQUILIBRIUM  AND  TRANSLATIONAL  EQUILIBRIUM 


We  now  consider  the  case  of  rotational  nonequilibrium  and  transla¬ 
tional  equilibrium  (Rj.  -»•  ®)  but  assume  R^.  »  1.  In  this  limit,  there  is  a 
single  longitudinal  mode  at  line  center  for  each  rotational  energy  level  J 
that  can  maintain  threshold  gain. 


The  assumption  R^.  “  implies  an  equilibrium  line  shape 

-X^  - 

^ 


(C-la) 


that  replaces  Eq .  (4a)  due  to  violation  of  Eq .  (3d).  Assuming  Av^  <<  Av^ , 
the  gain,  from  Eq.  (4c),  is 

-X^  - 

(2/71)  Gj(v)  =  e  Nj 

The  relation  between  threshold  gain  and  rotational  number  density  N^, 
for  those  rotational  energy  levels  that  can  support  a  single  longitudinal 
mode  at  line  center,  is 


(2/ti)G  =  N~ 

c  J 


(C-lb) 


For  the  case  of  laminar  mixing,  Eqs.  (4b)  and  (4d)  become,  respectively, 


Av  R 

'd  33-  -  ‘  -  2^*,  I  le  R-  L 

J  dr 


,  ,  ,,I/2  ‘“h  h 

itr— Ij-f 

dr  J,ra  N 


i-2L-  1  >  0 


(C-ld) 
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APPENDIX  C 

ROTATIONAL  NONEQUILIBRIUM  AND  TRANSLATIONAL  EQUILIBRIUM  (continued) 


where 

2/-  --  =  2D(2j)  -  (C-le) 

♦  r  =  (C-lf) 

Eq.  (C-lc)  can  be  integrated  upon  specification  of  z^,  6^,  and  R^.  The 
summation  in  Eq,  (C-lc)  is  evaluated  using  Eq.  (C-ld)  and  includes  those 
values  of  J  for  which  Eq.  (C-ld)  is  nonnegative.  Net  output  power  from 
each  rotational  trainsition  and  net  laser  power  are  found  from 


1/2  ^'’h  ^t 


^D^J,e  ^  ^r  J"  (’’^"2)  '  p"  Ij  zdz 


z . 
1 


d  r 


!„P  =  y  z_P , 

D  e  j  D  J,< 


(C-2a) 


(C-2b) 


Numerical  results  are  presented  in  Table  III  and  Fig.  9  for  the  case 


^i 


F  = 


R_  = 


0,001  -  0.5 

1/80 

0.0957 

100 


(J„  =  6) 


(^'j.m  = 


(C-3) 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation’s  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  Is  the  technical  staff’s  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory;  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermoraechanlcs,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exciraer  laser  development  Including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of -f ield-of “View  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors^  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  cooip'^und  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  ^olid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  mlcrowave/mllllmeter  wave  measurements, 
diagnostics  and  radloraeCry,  microwave/miHimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth’s  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


